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The color-electric flux tube of Abelian-projected (AP) SU(2) lattice gauge theory in the maximally 
Abelian gauge (MAG) is revisited. It is shown that the lattice Gribov copy efi'ect in the MAG is 
crucial for the monopole-related parts of the flux-tube profiles. Taking into account both the gauge 
fixing procedure and the effect of finite quark-antiquark distance properly, the scaling property of 
the flux-tube proflle is conflrmed. The quantitative relation between the measured AP flux tube and 
the flux-tube solution of the U(l) dual Abelian Higgs (DAH) model is also discussed. The fltting 
of the AP flux tube in terms of the DAH flux tube indicates that the vacuum can be classifled as 
weakly type-I dual superconductor. 



I. INTRODUCTION 

An intuitive explanation why a quark cannot be isolated as a free particle rests on the assumption that the QCD 
vacuum has the property of a dual superconductor Analogously to electrically-charged Cooper pairs being 

condensed in normal superconductors, magnetically-charged monopoles would be condensed in the QCD vacuum, and 
the dual analogue of the Meissner effect could be expected to occur. In the result, for example, the color-electric 
flux connecting a quark-antiquark {q-q) system would be squeezed into a quasi-one-dimensional flux tube |^ . 

This configuration provides a linearly rising potential between the quark and antiquark such that a quark cannot be 
separated infinitely from an antiquark spending a finite amount of energy. 

Remarkably, lattice QCD simulations with 't Hooft's Abelian projection typically in the maximally Abelian 
gauge (MAG), support this picture numerically. The distributions of the electric field and of the magnetic monopole 
currents, which can be identified after Abelian projection, have been measured in the presence of a static q-q pair 
(represented by a Wilson loop) within pure SU(2) 0,13 and SU(3) lattice gauge theories Q . It has been shown that 
the electric fiux is confined in a dual Abrikosov vortex due to a monopole current circulating around the q-q axis, 
signalling the dual Meissner effect. More quantitatively, the London penetration length of the electric field has been 
studied systematically within SU(2) lattice gauge theory j^. These authors compared the penetration length in MAG 
with a gauge invariant definition of the flux-tube profile. They came to the conclusion that the penetration length 
is gauge independent. A large-scale simulation on a 32^ lattice at a single value oi (3 — 2.5115 within SU(2) lattice 
gauge theory has been performed next |lo|. applying the fine-tuned gauge fixing algorithm (a mixture of overrelaxation 
(OR) algorithm and a realization of simulated annealing (SA)) [lll | in order to fight the lattice Gribov copy problem in 
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MAG and applying a noise reduction technique, the smearing of spacehke hnk variables before constructing a Wilson 
loop . Through this study, the dual Ampere law, a relation between the monopole current and the curl of electric 
field, has been confirmed with high accuracy. 

What is interesting and suggestive is that these numerical results gave strong hints towards the existence of a dual 
Abelian Higgs (DAH) model (we often call this the dual Ginzburg-Landau (DGL) model) as an effective model to 
deal with the QCD vacuum |l2l Il3l Il4| and with quark-induced hadronic excitations of the vacuum. In particular, 
the DAH model has a flux-tube solution corresponding to a q-q system. The mechanism of the flux-tube formation is 
nothing but the dual Meissner effect. The DAH model essentially contains three parameters, the dual gauge coupling 
(3g, the masses of the dual gauge boson tub and of the monopoles m^. The inverses of these masses are identified as 
the London penetration length and the Higgs coherence length, respectively. These lengths determine the width of 
the flux tube. The so-called Ginzburg-Landau parameter is defined as the ratio between the two masses, k = m^/mg, 
where the vacuum for k < 1 (> 1) is classified as type-I (type-H) superconductor. 

Determining the parameters of the DAH model, based on the comparison between the most elementary flux-tube 
profile measured in non- Abelian lattice gauge theory and the fiux-tube solution of the DAH model, is expected to be 
an important source of information on the QCD- vacuum structure itself, which would be helpful to learn about how 
the vacuum can function as a dual superconductor. Moreover, it might be useful in view of possible applications of the 
model in more complex physical situations (to describe baryons, for example in the SU(3) case). The first motivation 
has been, more or less, common to the above mentioned works and, in fact, such quantitative level of investigation has 
been ..tempted in Refs. RWU For example, the lattice data have been fitted by „me approximate, analytical 
fiux-tube solutions of the DAH model, and in this way the GL parameter have been estimated. The conclusions have 
varied, ranging from the vacuum belonging to the borderline between a type-I and type-II superconductor, with k ~ 1, 
as claimed in Refs. to a classification of the vacuum as a type-I superconductor, with k < 1, in Ref. 10]. The 

reanalysis of the profile data in Ref. 10] by fitting them with a numerical solution of the DAH model has supported 
the case of k ~ 1 Q . 

However, the systematic analysis of lattice fiux-tube data in MAG itself was incomplete such that doubts have 

remained whether the resulting parameters represent physical reality. At least, one has to check several basic properties 

of the lattice fiux-tube profile before one can seriously discuss the implications of the extracted DAH parameters. To 

mention the first, since one applies the MAG fixing, the lattice Gribov copy effect should be controlled, where the 
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OR-SA algorithm would be helpful for this purpose 11]. At second, one should check the scaling property of the 
profile with respect to changing the gauge coupling /3 of the lattice simulation. At third, one has to inquire how the 
profile behaves as a function of the q-q distance. One also should know how to compare the lattice fiux-tube profile 
with the flux-tube solution in the DAH model, if the above mentioned quantitative analysis is of interest. 

In our previous paper Q] , which was mainly devoted to the duality relating non- Abelian lattice theory on one hand 
with the DAH model on the other, we have carefully studied the flux-tube structure in the U(l) DAH model and have 
confronted it with some related data from our corresponding ongoing SU(2) lattice gauge measurements in the MAG. 
These studies, which will be discussed in the present paper in much more detail, have been done using a 32"* lattice 
with the OR-SA algorithm and with the smearing technique as in Ref. Then, based on the Hodge decomposition 
of the Abelian Wilson loop into the electric photon and the magnetic monopole parts, we have found there that the 
Abelian-projected lattice flux tube consists of two components, the Coulombic and the solcnoidal electric field, the 
latter being induced by the monopole currents circulating around the q-q axis. All this was in full analogy to the 
structure of the DAH flux-tube solution. We have also found that the Coulomb contribution cannot be neglected 
for any flux-tube length practically accessible in present-day lattice studies. 

In this paper we are going to present all our results concerning the flux-tube profile within AP-SU(2) lattice gauge 
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theory, obtained in the MAG, in a more complete way in order to meet the above requirements. The strategy of our 
study has been the foUowing one. Measurements have been performed using a 32'* lattice at various f3 values {f3—2.3, 
2.4, 2.5115, and 2.6). At first, we have investigated the lattice Gribov effect by comparing the profiles obtained from 
the widely used OR algorithm and from the OR-SA algorithm. In the context of the latter algorithm, we have also 
investigated the dependence on the number of gauge copies under investigation. At second, we have studied the scaling 
property of the flux-tube profile by comparing the profiles from various /? values, keeping the physical q-q distance 
approximately the same. The physical scales, the lattice spacing a{f3) for different values of /3, have been calculated 
through the measurement of the corresponding non-Abelian string tension. Throughout the profile measurements, 
smearing has always applied to the spatial link variables before constructing a Wilson loop. This procedure is meant 
to extract the profiles which effectively belong to the ground state of a flux tube; we have checked the (in)dependence 
of the flux-tube profile on the temporal extension of the Wilson loop. This effect has not received the due attention in 
the previous studies in Refs. QIESQ- This needs to be checked carefully when the Wilson loop is used to represent 
a static q-q source and the ground state is of interest. This procedure finally helps to reduce the noise. In a final step, 
we have assessed the DAH parameters by fitting the lattice data against the numerical DAH flux-tube solution. For 
this flt, we have not used the inflnitcly long flux-tube solution as it has been done in previous analyses Isi loi lioi llst . 
It should be noted that the use of the inflnitely long solution would be suitable only for that part of the electric fleld 
which is induced by the monopole part of the Wilson loop with sufficiently large temporal length . For our purpose 
to assess the DAH parameters, in this way we have taken into full account the finite q-q length effect in the fit. 

The paper is organized as follows. In section |nl we describe the procedures how to measure the Abelian-projected 
SU(2) flux-tube profile in MAG. Section UTTI presents the numerical results. In section Hvl we describe the results 
of fitting the lattice profiles by the DAH flux-tube solutions. Section is a summary and contains our conclu- 
sions. Preliminary results of the studies summarized in the present paper have been presented at the LATTICE2002 
conference ^] and in Ref. (l^ . 



II. NUMERICAL PROCEDURES 



In this section, we explain how to measure the profile of a color-electric flux tube on the lattice within the maximally 
Abelian gauge (MAG). We also develop the strategy to achieve a systematic, more detailed study of the flux-tube 
profile which takes into account the effect of the finite q-q length properly. We restrict the explanations of the methods 
to the case of SU(2) gauge theory. 

The numerical study of the flux-tube profile begins with the simulation of non-Abelian gauge fields. A thermalized 
ensemble of SU(2) gauge configurations {[/^i(to)} is generated by simulating the standard Wilson action 

Ssm){U] = /^E E (l - tr[C/^.(TO)]| (2.1) 

using the Monte Carlo heatbath method. Here U^,y{m) S SU{2) are plaquette variables constructed in terms of link 
variables U^{m) G SU{2) as 

U^^im) EE C/^(m)[/,(m + /i){7t(m + i))Ulim) . (2.2) 
The inverse coupling is given by /? = 
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A. Maximally Abelian Gauge fixing 

We put all the equilibrium configurations into the MAG. One exploits the gauge freedom of the SU(2) link variables 
with respect to gauge transformations g(m) 

Uf,{m)g''{m + fi) (2.3) 
in order to achieve a maximum of the following gauge functional 

^[^'] ^W^^' {^^U^.i^W H^)} , (2.4) 

where V is the number of sites in the lattice. The set of g{m) S SU{2)/U{1) for all site m represents the MAG fixing 
gauge transformation defined on m. 

For the numerical task to find the "optimal" g{m), in the past mostly an overrelaxation (OR) algorithm has been 
used. However, as it has been pointed out in the work in Ref. the OR algorithm is prone to fall into the nearest 
local maximum of Eq. (|2.4|l although the absolute maximum is of interest. This is due to the existence of many local 
maxima, which is known as the lattice Gribov copy problem. The only way known before to reduce the risk of being 
trapped in a wrong maximum is to explore many such local maxima by repeating the OR algorithm, starting each time 
from a new random gauge copy of the original Monte Carlo configuration; the ensemble of C/^(m) corresponding to 
only the highest of the achieved maxima was then understood as the gauge-fixed ensemble. It has then been proposed 
to use a "simulated annealing (SA)" algorithm with a (final) OR algorithm in order to prevent very poor maxima 
from entering the competition between gauge copies We may call this the OR-SA algorithm, which we mainly 
apply to fixing the MAG in the present paper. In the SA algorithm, the functional R[U^] is regarded as a spin action 



FiS) = R[U^] = ^ E {Sim)U^{m)S{m + fi)Ul{m)} , (2.5) 

where S{m) — g\m)T3g{m) corresponds to spin variables. The maximization of this functional is achieved by 
considering the statistical system given by the partition function 



Z = E^ exp 

{S(m)} 



(2.6) 



with decreasing the auxiliary temperature r — > 0. Practically, we first prepare a thermalized spin system at a certain 
high temperature, which is decreased gradually according to some annealing schedule until sufficiently low temperature 
is reached (r w 0). Then, final maximization is done by means of the OR algorithms. Notice that this algorithm only 
succeeds to escape from the worst local maxima, such that the above-mentioned inspection of many gauge copies per 
Monte Carlo configuration cannot be avoided. 

After the MAG fixing lOJ, the SU(2) link variables C/^(m) = Uff^im) are factorized into a diagonal (Abelian) 
link variable u^{m) e U{\) and the off-diagonal (charged matter field) parts c^{m), c*{m) e SU{2)/U{1) as follows 



The Abelian link variables Uf^{m) are explicitly written as 

u^(m) = e^^-(™) (0^(m)e [-TT,^)) . (2.8) 
The Abelian plaquette variables are constructed from the phase 9^{m) as 

e^,^{m) = e^,{m) + e^{m + (i)-e^,(m + v)-e^{m) (6*^^(771) G h47r,47r)) , (2.9) 
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which can be decomposed into a regular part 9pi,{m) G [— tt, tt) and a singular (magnetic Dirac string) part nf^^{m) = 
0,±1,±2 as follows 

6l^^(m) EE 0p^{m) + 2Tmp^{m) . (2.10) 

The field strength is defined by Opu{m) — Of^^i'm) ~ '2-'i^TLpy{m). Following DeGrand and Touissaint ,18.] , magnetic 
monopoles are extracted from the magnetic Dirac string sheets as their boundaries 

fc^(m) = -]^ep^pcrdunp„{m + ji) (£1234 = 1) , (2-11) 

where kp{rri) — 0, ±1, ±2 and rh denotes the dual site defined by to = m + (1 + 2 + 3 + 4)/2. Note that the monopole 
current satisfies a conservation law d'^kp{7fi) = formulated in terms of the backward derivative d'^. 



B. Correlation functions of Wilson loops involving local probes 



To find the flux-tube profile, one needs to measure the expectation value of a local probe 0{m) with an external 
source as {0{m))j (in our case j corresponds to an Abelian Wilson loop). Based on the path integral representation 
of {0{m))j, it can be rewritten as the ratio of {WAO{m))o and {Wa)o, where the subscript means the expectation 
value in the vacuum without such a source 

{ WAO{m) )o 
{Wa)o 



{0{m))j = j^^^ . (2.12) 



The Abelian Wilson loop is defined in terms of the Up{m) 

Up{ni) = exp [zJ20p{m)]. (2.13) 

The local operators 0{m) that we need to describe the structure of the flux tube are an electric field operator 

i9ii{m) = i{9n{m) — 2Tmi4(m)) , (2-14) 

and a monopole current operator 

2mki{m) , (2.15) 

where only spatial directions i = 1,2,3 are of interest. To avoid the contamination from higher states as much as 
possible, we have inserted the local probe 0{m) at t — tm = T/2 to minimize the effect from the boundary of the 
Wilson loop at t = and T. The local field operators are then evaluated over the whole x-y midplane erected in the 
center of the spatial extension of the Abelian Wilson loop {z = = R/2). In other words, the coordinates of the 
local operator are to = {x, y, Zm,tm) running over the midplane of the flux tube between a quark and an antiquark. 



C. Decomposition of the Abelian Wilson loop 



In order to see the composed structure of the flux-tube profiles, it is useful to apply the Hodge decomposition to 
the Abelian Wilson loop, which allows us to define the photon and monopole Wilson loops. We have shown in the 
previous work Q for the electric fleld profile that the photon Wilson loop induces exclusively the Coulombic electric 
field while the monopole Wilson loop creates the solenoidal electric field. At the same time, concerning the monopole 
current profile, the photon Wilson loop is not correlated with the monopole currents, while exclusively the monopole 
Wilson loop is responsible for the monopole current signal. 



We explain the decomposition using lattice differential form notations The Abehan Wilson loop in Eq. 
is written as Wa = exp[«(0, j)], where d{Ci) and j{Ci) are the Abelian link variables and the closed electric current. 
The Hodge decomposition of Abelian link variables leads to 

= A-^AO ^ A-\dS + Sd)0 

= A-^dS9 + A-^6e + 2TTA^^Sn, (2.16) 

where the second and third terms are identified as the photon link {6^'' = A~^5d) and monopole link (6'™° = 27rA~^(5n) 
variables, respectively. We do not need to fix the Abelian gauge in order to specify the first term in the second line, 
since it does not contribute to the Abelian Wilson loop due to Sj = 0: {A~^d69,j) = {dA~^59,j) = {A~^S9, Sj) = 0. 
Note that A is a lattice Laplacian and is its inverse, corresponding to the lattice Coulomb propagator. We have 
used d9 ~ 9 + 2TTn (see, Eq. (|2.1(J|l ) at the last equality. Inserting Eq. H2.16|l into the expression of Wa, the Abelian 
Wilson loop is decomposed as 

WaU] = exp[z(f?P\ j)] • exp[i(r-, j)] = Wph[j] ■ Wmo[j] ■ (2.17) 

We call Wph and Wmo the photon Wilson loop and the monopole Wilson loop. They are separately used to evaluate 
the photon and monopole parts of the profiles. 

D. Smearing of spatial links 

The shape of the fiux-tube profile induced by a Wilson loop depends on its size, Rx T, where R corresponds to 
the q-q distance and T the temporal extension. This means that the profile is infiuenced not only by the ground state 
but also by excited states, when a Wilson loop is used as an external source. If one is interested in the ground state, 
e.g., for the comparison with the flux-tube solution in the three-dimensional DAH model, in principle, one needs to 
know the profile at T — > oo. However, it is practically impossible to take this limit due to the finite lattice volume. 

Smearing is a useful technique to extract the profiles which belong to the ground state of a fiux tube effectively 
even with a finite T. We then see remarkable noise reduction when the size of the Wilson loop is large. The procedure 
is as follows. Regarding the fourth direction as the Euclidean time direction, we perform the following step several 
(Ns) times for the spacelike Abelian link variables: 

exp[i6'i(m)] 1-^ a exp [16*^(771)] -I- ^ exp i{9j{m) + 9i{m + j) — 9j{m + i)) , (2.18) 

where i,j = 1, 2, 3 and a is an appropriate smearing weight. 

To find an appropriate set of parameters (a, Ns), one needs to investigate the T-dependence of several quantities 
like the ground state overlap and the q-q potential. The emerging shape of the profile also should be checked for the 
effect of smearing. A numerical example at /3 = 2.5115 is shown in Appendix 1X1 We notice that this procedure seems 
to have practical limitations which become visible in the flux-tube proflle measurement. For the profile extracted with 
Wilson loops of size R <T, it works very well with large class of the parameter set (a, A^^). We could easily observe 
T-independence of the proflle within the numerical error. On the other hand, for the Wilson loops of size R > T, a,n 
extremely fine-tuned parameter set is required for smearing. However, we did not spend full effort to fix it. 

III. NUMERICAL RESULTS 



In this section, we present numerical results of the fiux-tube profiles measured over the x-y midplane of the q-q 
system (separated in z direction) using the Abelian, photon and monopole Wilson loops. We are going to clarify i) 
the lattice Gribov copy effect associated with the MAG fixing procedure and ii) the scaling property. 
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TABLE I: The non-Abelian string tensions and corresponding lattice spacings a(/3) for [3 = 2.3, 2.4, 2.5115, and 2.6 estimated 
by the relation a(/3) = ^faijo^^hyl with ^/d^^hys = 440 MeV. 



p 




a{P) [fm] 


2.3 


0.144(3) 


0.170(2) 


2.4 


0.0712(5) 


0.1197(4) 


2.5115 


0.0323(4) 


0.0806(5) 


2.6 


0.0186(2) 


0.0612(5) 



The numerical simulations were done at /3 = 2.3, 2.4, 2.5115, and (3 — 2.6. The lattice volume was always V = 32^. 
We have used 100 configurations for measurements. We have stored them after 3000 thermalization sweeps, and they 
were separated by 500 Monte Carlo updates. To study (i), we have generated several numbers of gauge copies from 
a given SU(2) configuration by random gauge transformations, each of which has undergone the OR-SA algorithm 
in the process of MAG fixing. The SA algorithm itself is applied with the temperature decreasing from t = 2.5 to 
T = 0.01. After that the OR algorithm is adopted with a certain convergence criterion. As the number of gauge 
copies, we have chosen iV^ = 5, 10 and 20, and have stored the configurations which provide the best maximal value 
of the gauge functional Eq. (|2.4(l within these Ng. We have also stored the same number of configurations (=100) 
from the OR algorithm in the MAG fixing with the same stopping criterion as in the OR-SA case, where always the 
first copy has been accepted {Ng = 1). To study (ii), we have used the configurations from the sample based on Ng = 
20 copies. The Abelian smearing parameters have been chosen Ng — 8 for a — 2.0. With this choice, the temporal 
length independence of the profiles induced by the Abelian Wilson loop is achieved within errors, at least for R < T. 
The same procedure was also applied to the spacelike photon and monopolc link variables before constructing each 
type of Wilson loop. 

A. Fixing the physical scale and choosing the flux-tube lengths 

The physical reference scale, the lattice spacing a(/3), has been determined from the non-Abelian string tensions 
by fixing ^Jophys = \/o'l/o. = 440 MeV. The non-Abelian string tension has been evaluated by measuring expectation 
values of non-Abelian Wilson loops with optimized non-Abelian smearing. The emerging potential has been fitted to 
match the form V{R) = C — A/ R + (Tl R- The resulting (dimensionless) lattice string tensions and the corresponding 
lattice spacings a{(3) in units of fm are shown in Table. 

To compare the profiles from various (3 values, it is quite important to put data into groups close to almost the same 

physical q-q distance because the finite length effect of the flux-tube system has to be studied simultaneously Q|. 

One might naively expect that the flux-tube profile has a good translational-invariant property along the q-q axis so 

that the difference in length does not matter when one follows the change in lattice scale. However, as shown in our 
I — I 

previous work 16], the finite length effect is not negligible as long as the photon part of the profile still contributes 
to the total profile. In Fig. ^ we then plot, for the four /3-values at our disposal, the physical length r = Ra{/3) in 
units of fm for various choices of the integer lattice flux-tube length R. This information is taken into account when 
we study the scaling property of the flux-tube profile. 

B. Assessment of the lattice Gribov problem 

We investigate how the flux-tube profile depends on the lattice Gribov copy effect due to the MAG fixing. As shown 
in Ref. [l]|. the density of monopole currents (in vacuum) is sensitive to the gauge fixing procedure; it decreases when 
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FIG. 1: Estimate of the physical q-q distance r in units of fm as a function of lattice distance R for various l3 = 2.3, 2.4, 2.5115, 
and 2.6 used in this study. 



larger R[U^] is achieved. Therefore, we expect that the monopole-related part of the profile crucially depends on the 
quality of the gauge fixing procedure. 

In Fig. [2 we show the electric field and monopole current profiles as observed over the flux-tube midplane at 
P = 2.3 for W{R, T) = VF(3, 3). Profiles, both with the use of the OR and the OR-SA algorithms, are presented. The 
dependence on the number of gauge copies under exploration is also investigated for the case of the latter algorithm. 
The upper, middle, and lower figures are the profiles from the Abelian, the monopole and the photon Wilson loops, 
respectively. 

We observe that the electric field and the monopole current profiles (except from the photon Wilson loop) are 
overestimated if the OR algorithm is applied. A possible explanation of this behavior is the following. The correlation 
between the monopole Wilson loop and the monopole currents are enhanced artificially due to denser monopole 
current system owing to the imperfect gauge fixing, which results in a larger contribution to the monopole current 
profile. Then, the strongly circulating monopole current around the q-q axis induces a strong solenoidal electric field. 
In this way, the electric field profile is also overestimated when the OR algorithm is adopted. It is interesting to note 
that since the photon Wilson loop is not correlated with the monopole currents, the corresponding electric field profile 
is insensitive to the Gribov copy problem. Finally, the impact of the Gribov copy problem on the monopole part is 
inherited also by the total flux-tube profile measured by the Abelian Wilson loop. Notice that the number of gauge 
copies to which the OR-SA algorithm is applied does not drastically change the profiles compared with the change 
from OR to OR-SA algorithm. This suggests that the tentative maxima successfully anticipated at the end of SA 
algorithm do not strongly differ in the monopole density. 

In Figs. 13 01 and [S] we show the same plot as in Fig. [3 for other (3 values, correspondingly choosing the Wilson 
loops: at /3 = 2.4 for 1^(4,4), at /3 = 2.5115 for W{&,Q) and at /3 = 2.6 for 1^(8,8). Here, the physical sizes of the 
respective Wilson loops are approximately the same (0.48 fm x 0.48 fm). We find that with increasing (3 (approaching 
to the continuum limit), the difference between the profiles from the OR algorithm and the OR-SA algorithm becomes 
clearer, i.e. the effect of the Gribov copy problem becomes more significant. Since the continuum limit is of interest, 



9 



0.30 
0.25, 
0.20 
0.15 
0.10 
0.05 
0.00 
-0.05 



-r 



o OR 

□ OR-SA(Ng=5) 

A OR-SA(Ng=10) 

V OR-SA(Ng=20) 



e- 



o OR 

□ OR-SA(Ng=5) 

A OR-SA(Ng=10) 

V OR-SA(Ng=20) 



0.14 
0.12' 
0.10 - 
0.08 
0.06 
0.04 
0.02 - 
0.00 
-0.02 



0.14 

0.12* 

0.10 

0.08 

0.06 - 

0.04 . 

0.02 - 

0.00 ■ 

-0.02 - 




o OR 

□ OR-SA(Ng=5) 

A OR-SA(Ng=10) 

V OR-SA(Ng=20) 





- B 
1 


mm 
1 


1 


1 





2 


3 


4 


5 6 



O OR 

□ OR-SA(Ng=5) 

A OR-SA(Ng=10) 

V OR-SA(Ng=20) 



e- 



o OR 

□ OR-SA(Ng=5) 

A OR-SA(Ng=10) 

V OR-SA(Ng=20) 




p/a 



p/a 



FIG. 2: The profiles of electric field (left) and monopole current (right) at /3 = 2.3 for the VF(3, 3) quark-antiquark system for 
gauge fixing according to OR and OR-SA algorithms, respectively, and the dependence on the number of gauge copies in the 
OR-SA case. Non-integer radiuses appear due to ofT-axis distances from the fiux-tube axis. Upper, middle, and lower figures 
refer to correlations with the full Abelian, the monopole and the photon Wilson loops. 



one needs to take care of this problem as already emphasized in Ref. [2] 
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FIG. 3: The same plot as in Fig. Hat (3 = 2.4 for 14^(4,4). 
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C. Does the flux-tube profile satisfy scaling ? 

We investigate the scaling property for groups of q-q distances according to Fig. using the best MAG-fixed 
configurations with Ng — 20. We choose three sets of physical distances: 

• at r - 0.48 fm: from /3 ==2.4 with W{R,T) = W^(4,4), /3 =2.5115 with W{6,6), and (3 =2.6 with W^(8,8), 

• at r - 0.80 fm: from (3 =2.5115 with W{10, 6) and /? =2.6 with W{13, 8). 

• at r - 0.84 fm: from f3 =2.3 with W{5,3) and /3 =2.4 with W{7,4). 

Here, the physical size of the temporal extension of the Wilson loop is taken approximately the same among all sets. 
This choice is made to normalize the systematic uncertainty for the flux-tube profile which might come from finite T 
even after smearing, especially for Wilson loops having R > T, which reflect contributions from excited states. We 
did not attempt to include profiles from such Wilson loops into the fit in Sec. IIVI The first set is shown in Fig. 
The other two sets are plotted together in Fig. [3 

We find that both the electric and monopole current profiles measured at different (3 values from the interval 2.3 
to 2.6 scale properly for each of the three groups of q-q distances. The remaining minor differences can be blamed to 
small differences in q-q distance and uncontrolled smearing effects. We also observe the following properties. Although 
the rotational invariance around q-q axis is poor for small /3, it is recovered with increasing /3. The electric field profile 
from the photon Wilson loop is very sensitive to the change of the q-q distance. Clearly, the shape of the electric field 
profiles from the photon Wilson loop in Figs. El and Hare different; this part of the electric field drastically decreases 
with increasing q-q distance. On the other hand, the electric field profile from the monopole Wilson loop remains 
almost the same with increasing q-q distance. The difference of the electric field profile coming from the full Abelian 
Wilson loop for different q-q distance can be explained by the change of the photon contribution. The large error 
of the monopole current profile from the Abelian Wilson loop in Fig. [7| is due to the large size of the Wilson loop, 
13 X 8, at /3 = 2.6. The statistics is not sufficient in this case. However, it is interesting to find that the decomposition 
of the Abelian Wilson loop into the photon and monopole parts helps to see a clear signal even with a number of 
configurations which is normally used for smaller Wilson loops. 
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IV. FITTING WITH THE U(l) DAH FLUX TUBE 

In this section, we discuss the quantitative relation between the extracted AP-SU(2) flux tube and the classical flux 
tube of the dual Abelian Higgs model (the DAH flux tube) through a fit of the former profile by the latter. In the 
fit, we take into account both the electric field and the monopole current profiles simultaneously. 

A. The dual lattice formulation of the DAH model 



The DAH flux-tube profile is calculated within the dual lattice formulation of the three dimensional DAH model in 
order to mimic eventual lattice discretization effects in the fit l2ll|. The lattice DAH action is 



'DAH 



\ y: F.Mr + ^ E - e'^'(™)$(^ + + (|ci,(™)p _ 

■i<j 2—1 



(4.1) 



where Fij is the dual field strength 

F.,j (m) = (to) + Bj (m + 1)- B, (m + ])- B^ (to) - 27rEy (m) . (4.2) 

Biim) and <I>(to) are the dual gauge field and the complex- valued scalar monopole field. The electric Dirac string 
in the dual field strength reflects the actual length of the flux tube. For instance, for a straight flux tube along 
the z direction, S12 — 1 for all plaquettes penetrated by this flux tube, otherwise = j^]. This action contains 
three parameters: the dual gauge coupling j3g — l/g"^, the dual gauge boson mass to^ = ^/2gv, and the monopole 
mass — 2%/Au. Here v corresponds to the monopole condensate, and A is the self-coupling of the monopole field. 
Writing the masses in terms of g, v and A is more familiar in the continuum form of the DAH model. Note that in the 
lattice formulation, all fields and parameters are dimensionless. The type of dual superconductivity is characterized 
by the Ginzburg-Landau parameter, k = m^/ms- In this definition, the cases k < 1 (> 1) are classified as type I 
(type II) vacuum. 

The fiux-tube solution is obtained by solving the field equations. The equation for Bi(rn) is given by gj'"^) — 
(3gXi{m) — 0. Similarly, the equations for the monopole field are q^R^,^^ = f3gm%X^{m) — and gj^Fff) — 
Pgm'^X^ (m) = 0. The superscripts R and / refer to the real and imaginary parts of the complex scalar monopole 
field. The explicit form of the field equations Xi{m), X^'{m), and X^{m) are given in Appendix IbI To solve 
the field equations numerically, we adopt a relaxation algorithm a la Newton and Raphson by taking into account 
the second derivative of the action with respect to each field (21|. We iterate this procedure until the conditions 
EmELi(^»("^))^ < 0.0001 and Em{(^^("^))^ + {^^{rn))^} < 0.0001 are satisfied and the change of the action 
for one iteration step AS'dah < 0.001. Within the possible range of the DAH parameters, we find that the solution 
is well- converged. 

The strategy of the fit is as follows. We fit the flux-tube profile induced by the Abelian Wilson loop with the 
DAH flux tube. Since we want to use the T-independent profile and it can only be T-independent if i? < T, we are 
restricted for this purpose (profile) to R < T. The electric field and monopole current profiles of the DAH flux tube 
is y/^SijkFjk and ^J^Ki (see. Appendix IB|| . which are regarded as representing On and 2'Kki of the AP-SU(2) field 
profile (see, Eqs. (|2.14() and (|2.15() '). The DAH field profiles are calculated with the same spatial volume as used 
in the SU(2) simulation, namely 32'^, imposing the same periodic boundary conditions for all three directions. The 
length of the DAH flux tube is taken equal to that of the AP-SU(2) flux tube to be fitted. We extract the profile all 
over the midplane cutting the flux tube between the quark and the antiquark. The DAH lattice spacing is assumed 
to be the same as a(/3) of the SU(2) lattice. Once the fit has found the optimal set of dimensionless DAH parameters. 
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the physical masses are fixed with the help of a(/3). To seek the set of parameters which provides minimum x^, we 
use the MINUIT code from the CERNLIB. 

After getting the set of the DAH parameter, we check whether this set can reproduce the composed internal 
structure of the electric field profile as a superposition of the Coulombic plus the solenoidal field by applying the 
Hodge decomposition as in Eq. (|2.16() to the dual field strength. Each field strength can be constructed from the 
DAH photon link (B^'' = 27rA-i(5S) and the DAH monopolc Unk = A-'^SF), respectively. Note that the field 

strength from the photon (monopole) links describe the Coulombic (solenoidal) electric field. 



B. Fitting results 



We fit the flux-tube profiles at /3 = 2.5115 from W{3,6), W{A,6), W{5,6), and W{6,6), and at /3 = 2.6 from 
M^(4,8), VF(5,8), W{6,8), PF(7,8), and W^(8,8). Here, the physical length of the temporal extension of the Wilson 
loop for these /3 values is approximately the same, 0.48 fm. We did not attempt to fit the profiles from the Wilson 
loops with R > T, since they still contain the contribution from excited states (T dependence) even after the smearing 
(see, Sec. IHDi and Appendix^. In the fit, we have taken into account the data from p/a > 2 for /3 — 2.5115 and 
from p/a > 3 for f3 = 2.6 to certain maximum radii which provide the positive expectation values for the field profiles. 
We have checked that the DAH parameters emerging from the fit are rather insensitive with respect to restricting the 
fit range (further increasing the minimal radius). 

In Table im we summarize the parameters obtained by the fit. In Fig.|Hl we show how the AP flux tube is described 
by the DAH one using the profiles from VK(6,6) aX j3 ~ 2.5115 and from 1^(8,8) &i j3 — 2.6. One can see that the 
profiles from the Abelian Wilson loop are reproduced. Remarkably, the resulting DAH parameters also reproduce the 
composed internal structure of the AP flux tube as well. In this sense, the fit which takes into account the finite q-q 
distance works very well. In Fig. El we plot the fitting parameters as a function of the physical q-q distance, where 
the scale of the masses is recovered by using the SU(2) lattice spacing a{(3). The maximum physical q-q distance is 
around 0.5 fm. We find that the (ig becomes large as increasing the q-q distance, r = Ra{0), while the masses of the 
dual gauge boson and of the monopole are rather stable. The constant fit of the masses using the stable data for 
r > 0.3 fm provides 

rriB = 1091(7) MeV, (4.3) 
= 953(20) MeV. (4.4) 

The GL parameter is then found to be 

<rri 

K= ^ = 0.87(2) < 1, (4.5) 

which means that the vacuum corresponds to weakly type I in terms of the classification of the dual superconductivity. 
However, we have noticed that the change of the dual gauge coupling Pg as a function of r indicates that the vacuum 
cannot completely be regarded as the classical one. In fact, if one defines an effective Abelian electric charge based 
on the Dirac quantization condition as Ccs = 47r/g = 471^//^; this coupling shows an anti-screening behavior; Ccff 
becomes large with increasing r. The constant behaviors of and indicate that various widths of the AP flux 
tube, which are defined by the inverse of these masses (the penetration depth rrig^ and the coherence length m^^), 
do not depend on r. This is established at least up to a q-q distance r = 0.5 fm. 
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TABLE II: The q-q distance dependence of the DAH parameter 





R 








xVdof 


fit range 


2.5115 


3 




n 4633('23'l 


349n('43~) 


135/67 


p/a> 2 


2.5115 
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n 071 1 (a) 


n 45951'! 7) 


n 373S('7~) 


99 9/81 


p/a> 2 


2.5115 


5 




4485('29'l 


n 4n9n('6~) 


77.4/81 


p/a> 2 


2.5115 


6 


0.0840(8) 


0.4504(21) 


0.4091(3) 


186/97 


p/a > 2 


2.6 


4 


0.0719(9) 


0.3372(72) 


0.2284(371) 


35.7/91 


p/a > 3 


2.6 
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0.3295(36) 


0.2884(29) 
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0.3004(5) 


37.3/91 


p/a > 3 
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0.0907(14) 


0.3360(10) 


0.3081(2) 


75.5/105 


p/a > 3 





FIG. 8: Some examples of the fitting at /3 = 2.5115 for R = 6 (upper row), and /3 = 2.6 for R = 8 (lower row). The solid line 
is the DAH flux-tube profile (obtained by the fit). The dotted and dashed lines correspond to its monopole and photon parts 
(as predicted using the fit parameters) . 
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V. SUMMARY AND CONCLUSIONS 

The main aim of this paper has been to present the flux-tube profile data within Abehan-projected (AP) SU(2) 
lattice gauge theory in the maximally Abelian gauge (MAG) in a quality and sufficiently detailed in order to warrant 
quantitative discussions from the point of view of the dual Abelian Higgs (DAH) model. We have mainly studied i) 
the lattice Gribov copy effect associated with the MAG fixing procedure and ii) the scaling property {(3 independence) 
of the flux-tube proflle using a large lattice volume, 32^. During these investigations, we have always paid special 
attention to the composed internal structure of the AP flux tube. We have also carefully monitored the effect of 
smearing. 

(i) We have found that the ffux-tube profile is very sensitive to lattice Gribov copy problem in the MAG, in particular, 
the monopole-related parts of the profile are strongly affected. The monopole current profile is overestimated if one 
uses a non-improved gauge fixing algorithm or just one single gauge copy. Since we do not know the real global 
maximum of R[U^] (see, Eq. H2.4|l ). we cannot insist that our result is the final one. However, we have obtained 
significantly corrected profile data by virtue of the over-relaxed simulated annealing algorithm converging within a 
moderate number of gauge copies. 

{ii) We have confirmed the scaling property of the ffux-tube profile, which has been achieved by using the well- 
gauge-fixed configurations and by taking into account the finite q-q distance effect properly. In fact, the flux-tube 
profile is strongly depending on the size {R and T) of the Abehan Wilson loops W{R,T). At finite q-q distance R, 
the photon part of the flux-tube profile is crucially contributing to the total Abelian electric field measured at any 
distance from the external charges. 

Finally, we have investigated the effective parameters of the U(l) dual Abelian Higgs (DAH) model by fitting the 
AP ffux tube with the DAH ffux-tube solution. We have adopted a new fitting strategy; in order to obtain the 
DAH flux-tube solution, we have defined the DAH model on the dual lattice and have solved the field equations 
numerically with the same size of the spatial lattice volume (32'^) and with the same periodic boundary condition as 
in the AP-SU(2) simulations. In particular, we have fully taken into account the finite q-q distance effect. As a result, 
we also could reproduce the composed internal structure of the AP ffux tube in terms of the DAH one. We have 
found that the dual gauge coupling depends on the q-q distance, which becomes large with increasing q-q distance. 
In this sense, the comparison between the AP ffux tube and the DAH flux tube considered at the classical level is not 
perfect. On the other hand, we have found that the masses of the dual gauge boson and the monopole are almost 
constant as a function of the q-q distance up to 0.5 fm, which take values of 1100 MeV and 950 MeV, respectively. 
The Ginzburg-Landau parameter takes a value slightly smaller than unity, indicating that the vacuum is classifled as 
weakly type-I dual superconductor. 

We should mention that, while we have chosen the DeGrand-Toussaint prescription for the electric flux and the 
monopole currents featuring in the description of flux tube in the Abelian projected gauge theory, this choice is not 
unique. We have made our choice because of the possibility to relate then the AP theory to the DGL description 
Haymaker et al. [2^ . on the other hand, have proposed an alternative definition to satisfy the the Maxwell equations 
even at finite lattice spacing. It would also be interesting to study the differences of the measured flux-tube proflle 
with their definition. 

Suggested by an effective bosonic string description Q > it is expected that the width of the ffux tube broadens with 
increasing q-q distance. If such an effect exists in the effective DAH description, too, we would have to see the change 
of the DAH mass parameters as a function of the q-q distance. Our results show that, at least until 0.5 fm, the width 
of the ffux tube is appearing as an almost stable vacuum property. It could be argued that the bosonic-string-like 
features of the ffux tube might become manifest only for much more elongated strings. In order to study the existence 
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of string roughening, one might be forced to study the profiles correlated with Wilson loops of much larger size. For 
larger R, if one takes small T, one can get the signal of the flux-tube profile with the help of smearing techniques. 
However, one has to care that such a profile does not immediately correspond to the physical profile at T — > oo. In 
order to check if the flux tube becomes broader, the profile should be investigated in a T-independent regime. 

It would be worthwhile to extend the strategy of this paper to the AP-SU(3) flux tube in order to discuss the 
quantitative relation between SU(3) gluodynamics and the U(l)xU(l) DAH model. 
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APPENDIX A: FIXING THE SMEARING PARAMETERS 



The smearing procedure for the spatial link variables as indicated in Eq. (|2.18|l successfully reduces the contribution 
from excited states. To find a set of optimized smearing parameters, the weight a and the number of smearing sweeps 
Ns, we need to investigate the behavior of the ratio 

Co = [Wa{R. T)f+'/[WA{R, T + 1)]^ (Al) 

as a function of R for some fixed values of T 1X\. For interpretation, we notice that this ratio turns into the ground 
state overlap in the limit T oo. We apply the smearing step repeatedly on the spatial links until we get a good 
ground state overlap. In Fig. 1101 we show the typical behavior of the ground state overlap before and after smearing 
dX (3 — 2.5115. We see clearly that Cq got enhanced to Cq ~ 1 as the result of smearing. This justifies to select, for 
this [3, the smearing parameters a = 2.0 and Ng — 8 which have led to the improved ground state overlap. Using this 
set, we can also confirm a typical improvement of the Abclian potential 

V{R, T) = HWa{R, T)/Wa{R, T + 1)) . (A2) 

This is shown in Fig. ^] where the potential is plotted as a function of T for some fixed R. After smearing, we see a 
clear pattern of plateaus ranging from T = 1 to large T, the height of which corresponds to the values of the potential 
V{R) at T oo. 

In Fig. ^1 we show the effect of the smearing to the behavior of the flux-tube profile at f3 = 2.5115 for 1^(6,4), 
W{6,6), and 14^(6,8). Before smearing, the shape of the electric field and monopole current profiles are depending 
on the size of the Wilson loop; the profile from the smaller Wilson loop is enhanced. After the smearing, we see a 
reduction of the electric field for all cases and at least, the profiles from W{6,6), and 1^(6,8) coincide within the 
numerical error (the error is also reduced). This suggests that T-independence of the flux-tube profile is now achieved. 
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FIG. 10: The ground state overlap before (left) and after (right) smearing at /3 = 2.5115. 




FIG. 11: The potential before (left) and after (right) smearing at /3 = 2.5115. 

The shape of the monopole current profile is insensitive to smearing. However, we find a remarkable reduction of the 
error. Note that the profile from W{6, 4) still did not converge into the same profile, which means that the chosen 
smearing parameter set is not adequate in this case. 

APPENDIX B: FIELD EQUATIONS OF THE LATTICE DAH MODEL 



We note the field equations of the lattice DAH model. For the dual gauge field, 

dS 



dBi{m) 
where 



Xi{m) = Fij{m) + Fji{m- j) + Fik{rn) + Fki{m- k) - Ki{m) 
The last term corresponds to the monopole current 



Ki{m) = -m% 



$^(m) (^^^{m + i) sin Bi{m) + ^\m + i) cosBi(m)^ 



-$^(m) (^^^{m + i) cos Bi{m) - + i) s\nBi{m)^ 



(Bl) 



(B2) 



(B3) 
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FIG. 12: The flux-tube profiles before (left) and after (right) Abelian smearing at /3 = 2.5115 for R = 6 (fixed). Three cases 
of the temporal length, T = 4, 6 and 8, are shown. 



For the monopole fields 

as 



/3gm|xf (m) = , 



a$«(m) 
^=/3,m|^/(m) = 0, 



(B4) 
(B5) 



where 



3 . 

6^^{m) -^U^^{m + i) COS Bi{m) - ^\m + i) sin Bi{m)^ 
+ ^$^(m - i) cos Bi{m - i) + $^(m - i) sin Bi{m - | 



(B6) 



X\m) = 6^'{m) -^U<P^{m + i) sin Bi{m) + ^^{m + i) cos Bi{m)j 

i=l 

+ ^$^(m - i){-sinBi{m - i)) + $^(m - i) cosBi{m - i)j | 



-ml^\m) ($^(m)2 + ^\mf - l) 



(B7) 
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